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1 INTRODUCTION 



ABSTRACT 

We consider constraints on the planetesimal population residing in the disks of AU 
Microscopii, Beta Pictoris and Fomalhaut taking into account their observed thick- 
nesses and normal disk opacities. We estimate that bodies of radius 5, 180 and 70 km 
are responsible for initiating the coUisional cascade accounting for the dust production 
for AU-Mic, Beta-Pic and Fomalhaut's disks, respectively, at break radii from the star 
where their surface brightness profiles change slope. Larger bodies, of radius 1000km 
and with surface density of order O.Olg cm~^, are required to explain the thickness of 
these disks assuming that they are heated by gravitational stirring. A comparison be- 
tween the densities of the two sizes suggests the size distribution in the largest bodies 
is flatter than that observed in the Kuiper belt. AU Mic's disk requires the shallowest 
size distribution for bodies with radius greater than 10km suggesting that the disk 
contains planetary embryos experiencing a stage of runaway growth. 

2 SCALING ACROSS THE COLLISIONAL 
CASCADE 



Recent visible band images taken with the Advanced Cam- 
era for Surveys on the Hubble Space Telescope well re- 
solve the vertical scale height of two edge on debris 
disks, the 12Myr old l|Barrado v Navascues et al] 1 19991 : 



IZuckerman et aij |200l'') dusty circumstellar disks of the 
MlVe star AU Micr oscopii (AU Mic) and the A5V star 
p Pi ctoris {(3 Pic) l|Krist et al.l l2005l : IColimowski et al.1 
l2006l ). Also resolved is the inner edge of Fomalhaut's ec- 
centric ring also allowing a measurement of the disk scale 
height (|Kalas et al.ll200a ). The vertical scale height, H, is 
related to the inclination dispersion of dust particles and so 
allows an estimate of the velocity dispersion of the smallest 
particles. The velocity dispersion of planetesimals sets the 
energy of inter-particle collisions and so affects a calculation 
of th e dust production rate through a collisional cascade 
(e.g. lKenvonll2002| : IWvatt fc Dentll2002l : iDominik fc DecinI 
l2003l : IWvatt et all I l2007l ). The velocity dispersion is also 
sensitive to the presence of larger bodies in the disk as 
gravitational scattering or stirring c auses an increase in 
the velocity dispersion wi th time (e.g.. IStewart fc Idall2000l : 
iKenvon fc Bromlevll200ll ). Here by combining observations 
of observed vertical thickness with estimates for the dust 
production and gravitational stirring rates we will place con- 
straints on the underlying planetesimal population in these 
disks. Because of the difficulty in resolving vertical struc- 
ture, previous cascade calculations have not used a velocity 
dispersion consistent with that estimated for these disks or 
estimated the role of gravitational stirring. 



We consider three disks with resolved vertical scale heights. 
The properties of these three systems along with the quanti- 
ties we estimate from them are listed in Table[T] For AU Mic 
and (3 Pic we list properties at the radius, r, from the star 
where there is break in the surface brightness profile. For 
Fomalhaut, we list properties in the ring edge. One of the 
observed quantities is the optical depth, f (A), at wavelength, 
A, normal to the disk plane. Because the absorption or emis- 
sivity coefficient of a dust grain with radius a is reduced for 
A > a, and there are more dust grains with smaller radii, we 
expect the optical depth to be related to the number density 
of particles of radius a ~ A (e.g., see discussion in section 4 
bv lWvatt fc Dent|[2003 '). As we only detect the dust parti- 
cles in scattered light or in thermal emission, we use scaling 
arguments to estimate the number of larger bodies residing 
in the disk. 

Another observed quantity is the disk thickness that we 
describe in terms of a scale height H that here is a half width. 
The disk aspect ratio is the scale height divided by radius; 
h = H/r. A population of low inclination orbits has {z^) ~ 



\/2h. Here i — 



^) and {i^) is the inclination 



dispersion. Subsequently we also refer to e 



where 



(e ) is the eccentricity dispersion. We assume a Rayleigh 
distribution of particle inclinations and eccentricities. 

We review how the dust opacity and disk thickness can 
be used to estimate the planetesimal size distribution. Dust 
production in a destructive collisional cascade can in its sim- 
plest form be studied with a power law size distribution. 



2 



The single power law form for the size distribution is in part 
based on the simplest assumption that the specific energy 
(kinetic energy per unit mass), QJj, required to catastroph- 
ically disrupt a body is a fixed number independent of body 
radius; (often 2 x l O^erg g"^ for icy bod ies is used based 
on the estimates bv lKenvon fc Luu1ll999l ). The number of 
particles with radius a in a logarithmic bin of size din a is 
predicted to be 

dN 



N{a) oc 



(1) 



din a 

Using a logarithmic bin gives the same scaling with 
a as a cumulative dist ribution Nya (see appendix A 
by lO'Brien fc Greenbergj SoOS). In an infinite destruc- 
tive self-similar collisi onal cascade, the exponent is pre- 
dicted to be g = 3.5 (ll)ohnanvil ll968l : iTanaka et all 1 19961 : 
iDavis fc Farinellalll997r iKenvonI I2OO2I ). The main asteroid 
belt, if fit wi th a single power law, has a lower exponent 
of q ~ 2.3 (jlvezic et aj l200lh . It is coUisionally evolved 
but deviates from q = 3.5 because of additional removal 
mechanisms (e.g., Yarkovsky drift and resonances) and be- 
cause the material propertie s depend non trivially on size 
ijO'Brien fc Greenberd l2005l ) . In contrast t he larger bodies 
in th e Kuiper belt are consistent with g ~ 5 (|Bernstein et al.l 
I2OO4I ). Because of their low number these do not collide 
often enough to be part of an ongoing destructive colli- 
sional cascade. The high exponent probably reflects condi- 
tions during the early solar system when planetesimals were 
growing as well as co lliding (e.g.. iWetherill fc Stewart. 1993i : 
iKokubo fc Mall 19961 ). 

The number of objects of radius a can be estimated 
from another of radius ad using the scaling relation 



N(a) = TV, 



\adJ 



(2) 



This relates the number of larger particles to the smallest 
and so observable particles. Estimates of the number of dust 
particles, Nd, as a function of their radius, ad, can be made 
from studies of optical, infrared and submillimeter observa- 
tions. It must be kept in mind that because of the uncer- 
tainty in the exponent q, it is difficult to be accurate when 
extrapolating over orders of magnitude in the size distribu- 
tion. 

The fractional area covered by particles of radius a or 
r(a) in a log radial bin can be similarly estimated. Because 
the opacity depends on the number per unit area times the 
cross section area, our assumed power law gives for the opac- 
ity integrated over a log radial bin 

dr 



r(a) 



Td 



\adJ 



(3) 



d In a \ad , 

where Td — na^siad) and s{ad) is the number of particles 
per unit area with radius ad in a log radial bin. Likewise the 
surface mass density 



E(a) = E 



^(-) 

\adJ 



(4) 



where ~ Tdpdad- For q = 3.5, most of the disk mass is 
in the largest particles or at the top of the cascade. Gravi- 
tational stirring and dynamical friction heating and cooling 
rates are proportional to the product of the surface density 
time the mass (e.g., equations 6.1 and 6.2 bv lStewart fc Idal 
I2OOOI ). scaling as 



E(a)m(a) = 



Edmd ( — ) 



7-9 



(5) 



where Edmd « TdPd^d- Even when the size distribution is as 
steep as that for the large objects in the Kuiper belt (q ~ 5) 
gravitational stirring is dominated by the largest bodies. 

The optical depth is related to the collision time. For a 
population of identical objects the collision timescale 

t,oi ~ {iTQ.)-\ (6) 

jHanninen fc Saloll 19921 ) where Q, is the mean motion (angu- 
lar rotation rate for a particle in a circular orbit) at radius 
r. Since the collision lifetime is proportional to the inverse 
of the optical depth, the timescale for a particle of radius a 
to hit another with the same size scale (again in log radial 
bins) is 



tcol,s{a) ~ tcol.d ( ) 

\adJ 



q-3 



(7) 



As ex plored by iDominik fc DecinI (120031 ): IWvatt et all 
(|2007| ). smaller particles are capable of dispersing a larger 
one if the specific energy of the collision exceeds the crit- 
ical value. The collision lifetime is shorter by a factor 
of ~ e^"*^ (Equa t ion 2 1,22 and associated discussion by 
iDominik fc DecinI |2003| ). where e a is the radius of a 
smaller particle capable of disrupting one with radius a. 
The parameter e is estimated by considering what energy 
projectile object can disrupt the target. 



7,2 



1/3 



(8) 



(approximating Equation 25 by IDominik fc DecinI l2003^ n. 
where w^^j is the relative or inter-particle velocity disper- 
sion. We expect the relative velocity dispersion is twice the 
particle velocity dispersion or v^^i ~ 2u^ . 

We can now estimate the collisional lifetime for par- 
ticles in a log radial bin taking into account collisions with 
smaller particles. After multiplying by Equation|51 Equation 
[7] becomes 



( ay I ^* 
\ad) \Q*D 



(9) 



consistent with 



For q — 3.5, the timescale tcoi(a) oc 
Equation 23 bv Dominik fc De cin (2003). The maximum ra- 
dius object that will disrupt during the lifetime of the system 
is found by setting tcoi{a) to the age of the system, tage, and 
solving Equation [9] fo r a. This estimate was also used by 
IWvatt fc Deii^ (120021 ) in their section 5.3. In other words 
we define a radius, atop such that tcoi{atop) ~ tage or 



atop — ad 



(tageSTdil)' 



For q — 3.5 this gives 

(^)'(6 



at 



ad 



u 

Oh 



(10) 



(11) 



^ The square root term in equation 25 by IDominik &: DecinI 
1I2OO3I) should be positive 
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where P is the rotation period at radius r. If the disk is hot- 
ter or older then a higher surface density disk that contains 
more massive bodies is required to initiate the colhsional 
cascade and account for the dust production. 

Objects of radius atop are those hkely to be currently 
initiating the colhsional cascade. Using Equation|4]with atop 
we can estimate the total surface density in these massive 
objects. As the disk grinds up and is depleted, more massive 
but lower number density objects can enter and generate the 
cascade. 



2.1 In relation to observables 

We first relate the disk aspect ratio, h, to the velocity dis- 
persion and the inclination and eccentricity dispersions. A 
population of low inclination orbits has {z^) ~ so 
i ~ \/2h. An is otropically scatteri ng disk is expected to have 
i ~ e/2 (e.g., Ilnaba et al.ll200ll ). At low eccentricity, the 
radial velocity dispersion is {v^) ~ (e^)«^/2, and the tan- 
gential and vertical velocity dispersions are iv^) ~ {v1) ~ 
(e^)t;|-/8, where vk is the velocity of a particle in a circu- 
lar or bit (e.g., see equations C10a,b bv lWetherill fc Stewart! 
1 19931 ). The total velocity dispersion is the sum of the three 



velocity components corresponding to u 



\{e>l or 



V3 



(12) 



These approximations are consistent with v^^, = (1. 2 5e -I- 
■^^^ vlr used by previous studies (jWvatt fc DentI l2002l : 



IWetheriU fc Stewartl[l993 ). 

In Equation |3] we described the scaling of opacity in a 
log radial bin. The normal disk opacity inferred from obser- 
vations at wavelength A, depends on the disk emissivity or 
absorption coefficient (here denoted Q) 



f(A) 



^Q{X,a)da. 



(13) 



This is consistent with our definition for T(a) (Equation [3]) 
and approximations com monly used in interpret ing observed 
fluxes (e.g., equation 1 bv lBackman et ahll 19921 relating dust 
opa city to flux and the definition given in the caption of Fig. 
6 bv lPantin et al.iri997l ) . The simplest models for the absorp- 
tion or emissivity coefficient of a particle estimate that these 
coefficients are 



Q{X,a) 



(14) 



1 for A ^ a 

(I)-" for X>a 

fe.g.. lB"ackman et al ] |l992l : IWvatt fc De"ntll2002l ) with n ~ 1. 
For n ~ 1 and q ~ —3.5 by integrating Equation [13] we find 
that r(A) ~ 4r{a = A). 

More detailed mode ling of the absorption coefficients 
(e.g., IPoUack ct all 1 1994 ) shows deviations from this sim- 
plest model with strong structure at specific wavelengths 
such as the lO/xm silicate feature. In addition, the expo- 
nent q, describing the dust size distribution, may not be 
well constrained or may not be the same for small dust 
particles as for larger ones (e.g., lAueereau fc Beusd [ioosl ) 
or for cm or m sized bodies. The wavelength at which the 
absorption coefficient begins to drop for equation [14] may 
de pend on dust compos ition (see discussion in appendix D 
bv lBackman et al.l[l993 ). Multi-wavelength observations are 



required to better model the size distribution and composi- 
tion of the dust. To take this uncertainty into account we 
describe our estimates in terms of a factor fr, such that 



-{a = A) = r(A)//, 



(15) 



that relates the opacity estimated at a wavelength based on 
observations to the size distribution of particles with radius 
equal to that wavelength. 

An estimate of normal disk opacity at a par- 
ticular radius requires modeling the surface brightness 
distribution IC olimows ki et al.l [2006'; 'K rist et aLF l2005l : 
lAueereau fc Beust 2006). Unfortunately, normal disk opac- 
ity estimates are available only at a few wavelengths for the 
three disks we are considering here and not all of these are 
based on multi- wavelength models. While optical and near- 
infrared wavelength observations tend to better resolve the 
disks, they may not accurately predict the mm size distri- 
bution (e.g., see the discussion comparing the optical and 
near-infrared opacities to that predicted from t he submil- 
limeter for AU Mic bv lAugereau fc Beus^|2006i ). We sum- 
marize the existing observed optical depth measurements for 
these three disks in Table[T]and in the associated table notes 
but note that there is uncertainty in the conversion factor 
fr between the measured optical depths and the opacity 
function that we have use here, r{ad), the optical depth in- 
tegrated in a log radial bin of size 1 for dust particles of 
size ad = A. As the opacity of smaller grains is sensitive 
to the removal process as well as collisions it is important 
to use observed opacity that is dominated by particles that 
are not affected by rad iative forces (e.g., see discussion by 
iDominik fc Decinl[20o3 ). 

We now convert Equation [11] into a form more easily 
computed from observables. The observables are the disk 
aspect ratio, h and the normal disk opacity t(A) at wave- 
length A. The size of the objects initiating the colhsional 
cascade when q = 3.5 



ato 



5.4km 



A 



W/im 

Qh 



Ah 



\ lOOAU 



2 X lO^erg g-i 



/ tags \ ' / h 
l^lO'^yrj V0.02 



10-2 



(16) 



Because we have scaled with the inclination or aspect ra- 
tio instead of the collision velocity the exponent of r and 
Mt differ from but ar e consistent with equation 36 by 
iDominik fc DecinI (l2003l). The relation also differs fror ti pre- 
vious work ( Wvatt fc DentI |2002| : [Pominik fc DecinI l2003l : 
IWvatt et al.l 20071 ) because we have based our estimate on 



a collision time scaled from the face on disk opacity at a 
particular radius rather than the total fraction of starlight 
re-emitted in the infrared. 

Inserting our value for the atop into equation [5] yields 
an estimate for the total disk density. 



S(atop) 



0.0018 g cm" 



Pd 



Ig cm~ 



2 X lowers 



ta 



lOOAU 



lO^yr 



h ' 
a02. 



4 



lOfim 



10-2 



We have assumed here that the coUision cascade started 
very early in the lile ol the system, however at early 
stages the inter-particle velocitie s were probably not high 
enough for destructive collisions l|Kenvon fc Bromlevll200ll : 
iDominik fc Decinlliooj ). If the timescale of the destructive 
cascade were smaller then atop and 'S{atop) would both be 
smaller than the estimates given above. 

The product of the density times the mass for the bodies 
initiating the cascade atop 



(Em) (atop) 



8.9 X lO^^g^ cm"^ 



(lOOAu) 



2 X 106 erg 
A 



I lovy 



(lO-O (4) 



lOfim 

(—V(- 



pd 



Lg cm~ 



(18) 



3 HEATING THE DISK WITH 
GRAVITATIONAL STIRRING 

We explore the idea that the observed thickness of the disk 
is due to gravitational stirring by bodies of mass, m^, surface 
density, Es, and size a^,. We define a mass ratio /i^ = 

and surface density ratio a a = ^ . If the disk is in coUi- 
sional equilibrium then we expect that e ~ 2i. 

In the dispersion dominated regime, and assuming that 
the dispersions of the tracer particles exceed those of the 
massive particles doing the stirring [i > is and e > es) 



1 d{i^) 
Q. dt 



s^isBJ,(l3)l3 



(19) 



(based on Equation 6.2 bv lStewart fc Idal l200d) w here /3 = 
I ~ 0.5 (corresponding to Equation 2.11 by Stewart fc Idal 
I2OOOI '). The function described by IStewart^Idal (|2000l ') 
Jz(/3 = 0.5) ^ 2.0. The coeffic ient B ~ 2InA and we es- 
timate A using Equation 2.7 bv lStewart fc Idal (|200G| ') 

3a*J'?. (20) 



A 



As the coefficient, B, only depends logarithmically on A we 
can use the scale height estimated from observations to esti- 
mate A and we can solve Equation [19] finding that i oc i"^''*, 
specifically 



i{t) 



2 In Anitas jJ.s 



1/4 



The above Equation can be inverted at time tage 



I In A^ta 



(21) 



(22) 



where we have set P to be the rotation period at r. In terms 
of observables this leads to a constraint on the largest bodies 



(17) 



with size ag 



2.4 X 10^*g^ cm ^ 



\0.02J \ 10^ 



\10Vr, 

m) (tooau) {12) ■ ^^^^ 



(Em)(a,) 



We note that the constraint on the product of the surface 
density times the mass of the largest bodies is independent of 
the disk opacity. In contrast the estimates for the top of the 
collisional cascade (size of object and density) are sensitive 
to the dust opacity. 

3.1 Connecting the size distributions 

Equation[Tn]gives us an estimate for the size of the bodies at 
the top of the coUisional cascade, and Equation [TS] gives us 
the surface density times mass in the disk for these bodies. 
This product is well below that needed to account for the 
disk thickness with gravitational stirring (Equation [23} . To 
find the size, as, of the bodies responsible for the gravita- 
tional stirring we must extend the size distribution beyond 
atop • 

Unfortunately, for bodies with sizes a > atop we can no 
longer assume a size distribution consistent with a collisional 
cascade. There are few guidelines on what type of power law 
to use for bodies greater than 10km. The only known sys- 
tem that differs significantly from the size distribution ex- 
pected from collisional evolution might be the largest bodies 
in the Kuiper belt that have size distribution with power law 
g ~ 5 l|Bernstein et al.ll2004f ). A variety of size distributions 
might be produced during the phase of planetesimal growth 
with low values for the expone nt q at the high mass end 



implying runaway grow th (e.g.. 



Kokubo & Ida 1996; Inaba et al 



WetheriU fc Stewart|[l993l: 



2OOII ). 



To place constraints on the size and density of the 
largest bodies and exponent of the size distribution for these 
bodies we compare our constraint on the product of the sur- 
face density and mass of the largest bodies to the surface 
density and size of the bodies initiating the collisional cas- 
cade. 

In Figure [1] we plot the constraint on the product of 
disk surface density times mass for AU Mic. This constraint 
corresponds to a surface density as a function of the radius 
of a body and is computed from Equation [23] using values 
listed in Table [T] and /r = 4. The horizontal axis is log 
radius instead of log mass so the slope of this constraint is 
-3. The conversion between mass and radius has been done 
with a density of 1 g cm~^. On this plot we have plotted as 
dotted lines two other constraints on bodies in the disk. We 
estimate that the most massive bodies cannot on average be 
closer together than their mutual Hill spheres. 



Efm) < 



' rriH 



(24) 



where the mutual Hill radius for two bodies of similar mass 
VmH = r (^g^) . This constraint gives the upper dotted 
line. We also require that the number of bodies not be ex- 
tremely low. 



E(m) > 



10m 



(25) 



This constraint is plotted as the lower dotted line. The range 



5 



E 




1 2 3 

logio a (km) 

Figure 1. The thick solid line shows the constraint on the prod- 
uct of the surface density times mass in the most massive bodies 
present for AU Mic, required to account for the disk thickness 
from heating by gravitational stirring. This is computed using 
Equation 1231 and values listed in Table [T] The upper dotted line 
shows the upper limit on the surface densities for these massive 
bodies set by requiring that they be on averaged spaced further 
apart than their mutual hill spheres (Equation 1241 . the lower 
dotted line shows the lower limit on their surface density set by 
requiring more than a few bodies of this mass reside in the disk 
(Equation |25}. The large circle is placed at the estimated location 
of the top of the coUisional cascade (computed using Equations 
1161 and 1181 and listed in Table 1). Arrows are shown with slopes 
predicted for size distributions with q = 3.5 and g = 5. The size 
distribution must connect the circle and the segment of the thick 
solid line that lies between the two dotted thin lines. 



of densities for the most massive bodies in the disk must lie 
on the solid one and between the two dotted ones. Also 
plotted on this plot is the estimated density, T,{atop), and 
radius, atop, of the particles initiating the cascade. Arrows 
are drawn for surface densities E(a) that have size distri- 
butions with exponents q = 3.5 and q — 5.0 and that have 
S(atop). The circle showing the top of the coUisional cascade 
must be connected to the thick solid line segment that lies 
between the two dotted ones to estimate the exponent of the 
size distribution for a > atop- 

The solid thick line segment between the two thin dot- 
ted lines in Figure [1] suggests that 1000km bodies reside in 
AU Mic's disk even though the coUisional cascade only re- 
quires bodies of radius a few km. We have checked that our 
estimated value of 12 for log A is consistent with the mass 
of these 1000 km bodies and the disk thicknesses, (equation 
I20p . For q > 4 most of the disk mass resides in the most 
massive bodies. Connecting the circle with the line segment 
requires a slope shallower than q = 3.0. Most of the disk 
mass must reside in 1000 km embryos in AU Mic's disk to 
account for its thickness even though only km sized bodies 
are required to account for its dust production. 

Figure [5] and O are similar to Figure [T] except computed 
for j3 Pic's and Fomalhaut's disks also using parameters 
listed in Table [1] We attribute the differences in these figures 
primarily to the observed thickness as atop oc h^'^^^ (equa- 
tion ll6|) . /3-Pic's disk is quite a bit thicker than Fomalhaut's 
or AU Mic's so its coUisional cascade is more efficient and 
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Figure 2. Similar to Figure [T] except for /3 Pic's disk. 
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Figure 3. Similar to Figure [T] except for Fomalhaut's disk. 



so requires higher mass progenitors. Fomalhaut is older al- 
lowing a lower density disk to account for the thickness. 

Gravitational stirring requires similar sized embryos for 
the three disks but for Fomalhaut the mass and surface den- 
sity of the bodies is only an order of magnitude larger than 
that predicted from estimating the top of the coUisional cas- 
cade. Nevertheless, the bodies we infer at the top of the 
coUisional cascade are not sufficiently dense and massive to 
account for the thickness of this disk. 

A comparison between the surface densities in the bod- 
ies required to account for the disk thickness and that pre- 
dicted at the top of the collision cascade allows exponents 
g < 3, 3.5 and 4.5 for the three disks AU Mic, /3-Pic and Fo- 
malhaut, respectively. The extremely shallow exponent for 
AU Mic at the top end suggests that the size distribution 
deviates from power law form. A curve in the size distribu- 
tion at the high mass end has been predicted by models and 
simulations of planetesimal accretion when the disk contain 
embryos in a stage of runaway growth ( Wetherill fc Stewart! 
1 19931 : iKokubo fc ldall 19961 : llnaba et al.ll200ll ). 

We have only considered the effect of gravitational stir- 
ring in the dispersion dominated regime. Now that we have 
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an estimate for the masses of the most massive bodies resid- 
ing in these disks, we check this assumption. Only for a > 
1.2 X lO^km does a body's Hill radius approach a scale height 
ri for an incUnation i = 0.01. The dispersion dominated 
gravitational stirring estimate used in Equation [19] (rather 
than a sheer dominated one) is therefor reasonable. Previous 
work has found that passage thr ough the sheer dominate d 
regime is comparatively fast (e.e:.. lKenvon fc Bromlevll200lh . 
A better estimate would take into account both regimes, 
though the improved constraints on the massive bodies 
should not significantly deviate from those estimated here. 



4 DISCUSSION 



We have used estimates of collisional cascades (e.g.. iKenvonI 



12002 ; Dominik fc Decin 2003; W vatt et al. 2007) to estimate 
the size and surface density of the bodies responsible for ini- 
tiating the collisional cascade. We have done this for 3 debris 
disks, that of AU Mic, /3-Pic and Fomalhaut, with resolved 
vertical structure estimating that these bodies have radii of 
4, 180, and 70 km, respectively. We have estimated these 
at the radius at which the surface brightness profile changes 
slope (also called the break radius). The bod y sizes are a few 
times larger than previous estimates (e.g., IWvatt fc DentI 
|2002| ). The differences arise because we have based our esti- 
mate on a collision time scaled from the face on disk opac- 
ity at a particular radius rather than the total fraction of 
starlight re-emitted in the infrared and we have used the ob- 
served disk aspect ratio to estimate the velocity of collisions. 

Assuming that the smallest particles are heated solely 
by gravitational stirring from the largest ones, the disk thick- 
ness can be used to place a constraint on the product of 
the surface density times mass of the largest bodies (Equa- 
tion [T8|. From this we infer that 1000km radius bodies or 
planetary embryos are likely to reside in these three disks. 
The large body sizes do not conflict with t he lack of ob- 
served gaps in the disks (|Quillenll200"6l . l2007l ) except possi- 
bly for the extreme high mass end allowed for /3-Pic's disk. 
A comparison between the surface densities in these bodies 
and that predicted at the top of the collision cascade al- 
lows exponents q < 3, 3.5, 4.5 for the three disks AU Mic, 
/3-Pic and Fomalhaut, respectively. The shallow exponent 
for AU Mic at the top end suggests that this disk con- 
tains embry os in a stage of runaway growth, as predicted by 
simulations (IWetherill fc Stewartll 19931 ; iKokubo fc Idalll996l ; 
llnaba et al.l 12001 ). For all three disks we infer that most of 
the disk mass is likely to reside in embryos and estimate that 
the surface densities are of order 10~^g cm~^. 

A number of simplifying assumptions went into estimat- 
ing the properties of the top of the cascade. We assumed only 
a single power law form for the size distribution, however, 
the specifi c energy for dispersion is predicted to depend on 
body size (|Benz fc Asphau^ll999l ) so a single power law is 
probably not a good assumption. The disks may not have 
been sufficiently excited for efficient dus t production dur- 
ing t he entire lifetime of these systems (jPominik fc DecinI 
l2003l ). A shorter collisional lifetime would lead to a lower 
surface density and size estimated for the top of the cascade 
(see Equations 1171 116|l , though taking into account the de- 
pendence of the specific energy on size in the regime where 
self-gravity is important would increase the surface density 



of larger bodies and might decrease the size at the top of 
the cascade. The sizes at the top of the cascade predicted 
here are nearing the threshold for a destructive equal-mass 
collision at a velocity estimated from the disk thickness, par- 
ticularly in the case of Fomalhaut that has a very thin disk 
but has a large estimated atop- 

Our estimate of the gravitational stirring rate neglected 
the role of dynamical friction from smaller particles and the 
sheer dominated regime. Both should be taken into account 
to improve the estimate of size and number of the largest 
bodies residing in these disks. 

Better modeling of the dust distribution using multi- 
wavelength observations and high angular resolution imag- 
ing would significantly improve constraints on the small ra- 
dius end of the size distribution. While we have found nor- 
mal disk opacity measurements in a few wavelengths in the 
literature, the different wavelength estimates, different as- 
sumptions for the assumed size distributions and different 
procedures for modeling the data make it difficult to con- 
strain and compare the dust size distributions and normal 
disk opacities among the disks. 

We have discussed ways to improve the estimates in- 
troduced here. We now discuss possible implications based 
on these predictions. If the size distributions inferred here 
are common then longer lifetimes would be predicted for 
dust production because the larger bodies (inferred here), 
entering the cascade later, contain a reservoir of mass avail- 
able for dust production at later times. The distribution of 
disk properties as a function of age can be used to place 
constraints on planetesimal growth models as well as dust 
production. 

We have only considered opacities at particular radii 
for these disks. For AU Mic and /3-Pic we chose radii at 
which there is a break (or change in slope) in the surface 
brightness profile. If the disk aspect ratios do not strongly 
vary with radius then Equation [23] implies that the product 
of the mass times the surface density in the largest bod- 
ies, T,m{as) oc r~^^^ is only weakly decaying with radius. 



Compare this to T,{ato 



-'1^ and atop oc r~V-"/3 
predicted via Equations I16I17I Both Yj(atop) and atop must 
drop rapidly with radius. If disks are not extremely thin 
at larger radii then either there is another source of heat- 
ing at large radii accounting for the disk thickness, or dust 
particles detected at large radii originate from inner radii 
and are either blown out or are on highly eccentr ic orbits 
(|Augereau fc Beustll200^ ; IStrubbe fc Chiangll2006D . A thin 
and sparse disk will not efficiently produce dust as the col- 
lisions are not destructive. Consequently multi-wavelength 
observations resolving disks as a function of radius should 
be able to test the utility of the estimates explored here as 
well as better probe planetesimal growth and evolution with 
radius. 
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Table 1. Debris Disks with measured thicknesses 



Stellar and Disk Properties 



Row 




AU Mic 


13 Pic 


Fomalhaut 


1 


M.(Mq) 


0.59 


1.75 


2.0 


2 


Age (Myr) 


12 


12 


200 


3 


r(AU) 


30 


100 


133 


4 


h 


0.019 


0.05 


0.013 


5 


T(A,r) 


3 X 10-3 


5 X 10-3 


1.6 X 10-3 


6 


A (/tm) 


1 


10 


24 



Estimated Planctesimal Properties 

7 atop (km) 

8 S(atop) (g cm-2) 

9 'Em{atop) (g^ cm-^) 

10 Em(as) (g^ cm-2) 



4 


180 


68 


0.00005 


0.005 


0.002 


1014.5 


1021.0 


lOis.s 


1024.1 


1026.2 


1022.7 



By Row. 1) Refer ences for the stellar masses: iHoudebine fc D oyle (1994'); 'Crifo ct al.' 1*1997'): 
Son^ et al.|| 200l|). respectiv ely. 2) References for the ages:|Barrado y Navascucs ct al. (1999); 
Barrado v NavascuesI lll998l) . 3) The radii are chosen to be whore there is a break in the surface 



bright ness profile as described by Kri st et al., (i2005i) ; Golimowsk i et all (|2006|) ; [Kalas et al.l 
^2003), respectively. 4) The aspect ratio h = H/r for H the half width half max of the disk 
at radius r. Aspect ratios are taken from the same references as the break radii listed in 
row 3. 5,6) The normal disk opacity f at wavelength A is given. References for normal disk 
opac ities: The normal d i sk opa city for AU Mic is estimated for Ifi sized particles from Fig. 
6 by lAugereau fc Beusti l|2006h based on im ages i n the optical and near infrared. That for 
/3-Pic is taken from Fig. 6 by |Pantinet al.' (1997) based on mid-infrared spectra. That for 
Fomalhaut is from Table 1 bv iMarsli^ t al. (2005) predicted for a reference wavelength of 
24fiiai based on 350, 160 and TOfira imaging. 7) The radius of objects initiating the coUisional 
cascade, atop is estimated using Equation 1161 8) The surface density Jl(atop) is estimated 
using Equation 1171 9) The product of the surface density times the mass {T:m){atop) is 
estimated for bodies initiating the coUisional cascade. 10) The product of the surface density 
times the mass is estimated using Equation [18] for bodies responsible for thickening the disk. 
Computed quantities listed in rows 7-10 have been done with parameter /t = 4 (defined in 
Equation [TSj . 



